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Dwarf nova oscillations and quasi-periodic oscillations in 
cataclysmic variables — VI. Spin rates, propellering, and 
coherence 



o 
o 
o 

(N 

> 

Q 

o 

6 



> 

in 
o 



o 



X 



Brian Warner^'^* and Magaretha L. Pretorius,^* 

^Department of Astronomy, University of Cape Town, Rondebosch 7700, Cape Town, South Africa 

^School of Physics and Astronomy, University of Southampton, Highfield, Southampton S017 IBJ, United Kingdom 



Accepted . Received ; in original form 



ABSTRACT 

We examine published observations of dwarf nova oscillations (DNOs) on the rise and 
decline of outbursts and show that their rates of change are in reasonable agreement 
with those predicted from the magnetic accretion model. We find evidence for propel- 
lering in the late stages of outburst of several dwarf novae, as shown by reductions in 
EUVE fluxes and from rapid increases of the DNO periods. Reanalysis of DNOs ob- 
served in TY PsA, which had particularly large amplitudes, shows that the apparent 
loss of coherence during late decline is better described as a regular switching between 
two nearby periods. It is partly this and the rapid deceleration in some systems that 
make the DNOs harder to detect. 

We suggest that the 28.95 s periodicity in WZ Sge, which has long been a puzzle, 
is caused by heated regions in the disc, just beyond the corotation radius, which are a 
consequence of magnetic coupling between the primary and gas in the accretion disc. 
This leads to a possible new interpretation of the 'longer period DNOs' (IpDNOs) 
commonly observed in dwarf novae and nova-like variables. 
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1 INTRODUCTION 

Dwarf nova oscillations (DNOs) were first recognized in out- 
bursts of dwarf novae and in nova-like variables; that is, 
during stages of high mass transfer rates (M) from com- 
panion to the white dwarf primary of cataclysmic variables 
(Warner & Robinson 1972). Their short time scales (~ 5- 
50 s) indicate originations close to the surface of the pri- 
mary, but their large variations in period distinguish them 
from the short period highly stable magnetic accretors such 
as DQ Her and AE Aqr (for a general review see Warner 
1995a). Paczynski (1978) suggested that these variations re- 
sult from the torque acting on an accreted equatorial belt on 
the primary during phases of changing M . This model has 
been extended in Warner (1995b) and the earlier papers in 
this series (Woudt & Warner 2002, Warner & Woudt 2002, 
2006; Warner, Woudt & Pretorius 2003; Pretorius, Warner 
& Woudt 2006); a review is given in Warner (2004). 
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In this contribution we reinvestigate some early observa- 
tions, in the light of later developments, looking specifically 
for additional evidence relevant to the Low Inertia Magnetic 
Accretor (LIMA) model (Warner & Woudt 2002; hereafter 
WW2002). In Section 2 we compare predicted and observed 
rates of change of spin on the rise and fall of dwarf nova 
outbursts; in Section 3 we look at possible incidences of a 
propellering phase late in dwarf nova outbursts; in Section 4 
we consider the question of 'coherence' of DNOs; in Section 

5 we suggest a model for the multiple DNOs in WZ Sge, their 
behaviour during outburst, and their possible relevance to 
the distinct group of longer period DNOs (IpDNOs). Section 

6 contains a short general conclusion. 



2 SPIN-UP AND SPIN-DOWN 

Our concept of the progress of DNOs through a dwarf nova 
outburst is that, at the beginning of outburst, material mov- 
ing inwards and reaching the inner parts of the accretion disc 
is either magnetically channeled into the equatorial regions 
of the primary or reaches the equator without channeling. 
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depending on the magnetic field in the inner regions. Even 
if the field of the primary is quite small, differential rota- 
tion of the ensuing equatorial belt can use this seed field 
and may amplify it to the point where channeling can occur 
(WW2002). The accretion zones generated by chaimeling 
are directly or indirectly responsible for the DNOs. Below 
some critical field of the primary (estimated in WW2002 to 
be ~ 10"' G) there will be no channeling and no DNOs dur- 
ing outburst. This is why otherwise very similar systems can 
exhibit completely different behaviours of their DNOs (e.g. 
VW Hyi and EK TrA; Warner & Woudt 2006). 

Moreover, for primary fields less than ~ 10^ G the outer 
regions are not coupled to the interior of the primary and 
so the accreted gas is free to rotate as an equatorial belt 
as it builds up in mass during the outburst. It is, however, 
magnetically coupled to the disc in the same way as in inter- 
mediate polars (IPs) , and therefore as the inner edge of the 
disc moves inwards through the increased thrust of the mass 
inflow, accessing shorter periods of rotation, the belt is spun 
up, and vice versa after maximum of outburst. Through the 
maximum phases of outburst the accretion is neax equilib- 
rium disc accretion, as is evidenced by the higher stability 
(e.g. observed low P) of the DNOs during those times. 

We here extend the quantitative aspects of this model, 
as developed in WW2002: we note that from the standard 
theory of accretion from a disc into a magnetosphere (e.g. 
Ghosh & Lamb 1979; Chapter 7 of Warner 1995a) the an- 
gular acceleration generated by the combined material and 
magnetic torques acting on an equatorial belt with mass 
M(6), given in equation (5) of WW2002, can be written in 
terms of the period of rotation P of the belt as 

P = 2.67 X lQ-^n{u)s)cjl^^MwM22\b)Mt^^{l)P'^/^ (1) 

where n{us) is the torque function, ujs is the "fastness pa- 
rameter" , given by where ro is the radius of the 
inner edge of the disc and Vco is the corotation radius, found 
from rj?o = GM{l)/n^{b), M{h) is the mass of the bch, Q.{h) 
is its angular velocity, M is the rate of mass transfer onto the 
primary, and M(l) is the mass of the primary. Subscripts 
denote cgs values in powers of ten, with the exception that 
Ml is used for meiss in solar units. 

An outburst can be partitioned into two sets of condi- 
tions: on the rise and on the fall the accretion process is not 
close to equilibrium, i.e. the rotation periods of the belt and 
the inner edge of the accretion disc arc significantly different, 
but around maximum the accretion is close to equilibrium, 
with the result that n(ws) ~ 0, ojs ~ 0.975 and P 0. On 
the rise and fall, however, we have n{ijJa) ~ 0.7 (Wang 1987) 
and 0.25 < ^ 1-0, so the factor ujV^ in equation (1) is 

1/3 

bounded by 0.63 < los ^ 1-0; the upper bound comes from 
the requirement that gas threading onto field lines must lie 
within the corotation radius, and the lower bound arises 
from (a) the requirement that ro ^ R(X) ^'id (b) the ob- 
served low values of P for the DNOs that imply that the 
inner edge of the disc is never very far from the surface of 
the primary. 

Next we note that M during the outburst phases of 
interest hero probably does not vary greatly. The models 
of accretion during outburst produced by Cannizzo (1993) 
show that, from the mid-point of optical rise to the midpoint 
of optical decline, M onto the primary varies only about a 



factor of two about the average. These correspond roughly 
to the phases over which most DNOs are observed. There- 
fore the mass of the equatorial belt (which we assume to be 
very small at the outset) during the interesting parts of the 
outburst can be approximated by M(6) oc T, where T is the 
time since the start of outburst. 

For the two non-equilibrium phases we therefore have 
(with appropriate choices of sign) 

P ~ ±1.6 X 10"* [Mi(l)/0.6]^/^ (P/25 s)^/^ T"' (2) 

where T is measured in days. 

We therefore expect P ~ 1 x 10~* on the rise and 
P ~ 3 X lO^'"" on the fall of a normal dwarf nova outburst. 
In Tables 1 and 2 we tabulate values of P deduced from 
published observations. There are relatively few usable ob- 
servations: for the rise this is because very few dwarf novae 
have been caught early on the rise, and for the decline it is 
because the DNOs often become too 'incoherent' to detect 
with Fourier transforms late in outburst. As an example, 
VW Hyi is one of the most intensively studied systems (see 
list of observations in Woudt and Warner 2002) but has 
never been caught showing DNOs on the rise, and there is 
only one detection of a DNO on the decline before the rapid 
deceleration phase begins, as described below. And we need 
to avoid the late rise or early fall where the low P character- 
istic of maximum still has infiuenco. However, it is immedi- 
ately obvious from Tables 1 and 2 that |P| is roughly a factor 
of five or more larger on the rise than in mid-decline. In the 
LIMA model this is expected because of the smaller iner- 
tia of the equatorial belt in the rising phase, compared with 
that of the almost fully accreted mass late in decline. Fur- 
thermore, after allowance for the range of parameters (e.g. 
Mi(l) ~ 1.2 and P ~ 10 s in SS Cyg, which have opposite 
infiuences in equation (2) - and this is a general rule, the 
higher mass primaries will generate shorter period DNOs) 
there is moderately good agreement with the predictions of 
equation (2). This is new evidence for the viability of the 
LIMA model. 

The notably larger values of P in Table 2, observed for 
OY Car, SS Cyg and TY PsA late in decline, are discussed 
in the next Section, along with a similar effect seen in VW 
Hyi. 



3 PROPELLERING 

As the system passes through the later declining phases of 
the outburst, M is rapidly reduced, acting to increase the 
radius ro of the magnetosphere, but at the same time any 
sheer-induced magnetic field in the belt is declining, which 
acts to reduce ro. The rate of increase or decrease of ro there- 
fore depends on details of these changes, which are currently 
unknown. One possibility is that ro increases much faster 
than Tco (which is determined by the spin period P{b)), 
catches up with it and turns an lJs ^ 1 accretor into an 
^ 1 non-accretor. That is, gas will be centrifuged out- 
wards on the field lines, in a propeller action, as has been 
suggested for VW Hyi (WW2002). 

In VW Hyi, unlike other dwarf novae, the late decline 
shows DNOs of large amplitude. These have a rate of in- 
crease of period nearly two orders of magnitude larger than 
predicted by equation (2), showing that some additional 
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Table 1. DNOs on early rise. The fifth column lists the range in time over which the 
period evolution given in the fourth column occurred. 



Star Outburst T (d) Pdno (s) Range (d) P/IQ-^ References 



SS Cyg 


Jun 1984 


1.9 


10.1 - 


9.2 


0.19 


8 


1 


SS Cyg 


Mar 1997 


1 


7.81 - 


6.59 


0.20 


7.1 


2 


AH Her 


Jun 1976 


2 


31.3 - 


30.7 


0.1 


7 


3 


AH Her 


Jun 1978 


2-3 


27.8 - 


25.0 


1.0 


3.2 


3 


SY Cnc 


Jan 1975 


2-3 


29.50 


- 26.71 


1.0 


3.2 


3 



References; 1. Jones & Watson 1992; 2. Mauche & Robinson 2001; 3. Patterson 1981. 



Table 2. DNOs during mid- and late-decline. The fifth column lists the range in time over 
which the period evolution given in the fourth column occurred. 



Star Outburst T (d) Pdno (s) Range (d) P/IO"^ References 



AH Her 


Jun 1978 


7-10 


24.25 


- 25.75 


3.0 


0.58 


1 


SY Cnc 


Jan 1975 


8 


27.80 


- 29.08 


1.0 


1.5 


1 






9 


29.08 


- 31.20 


1.0 


2.5 


1 


OY Car 


Nov 1980* 


12: 


19.44 


- 20.48 


1.0 


1.2 


2 






13: 


2048 


- 27.99 


1.0 


8.7 


2 


TY PsA 


Jun 1984* 


10 - 13 


25.2 - 


30.2 


3.0 


1.9 


3 






13 - U 


30.2 - 


- 38.0 


1.0 


9.0 


3 






U - 15 


38.0 - 


- 47.0 


1.0 


10.4 


3 




Sep 1984 


1-2 


26.5 - 


- 36: 


1.0 


11: 


3 


V1159 Ori 


Dec 1993* 


10 - 12 


31.61 


- 32.21 


2.0 


0.35 


4 






12 - 15 


32.21 


- 33.92 


3.0 


0.66 


4 


SS Cyg 


Sep 1984 


4 


7.8 






0.46 


5 






6 


10.4 






0.66 


5 






8 or 9 


10.7 




~ 0.5 


6-10 


5 



Notes: * indicates a superoutburst. Italicized entries axe for late-decline. 

References: 1. Patterson 1981; 2. Schoembs 1986; 3. Warner, O'Donoghue & Wargau 1989, 

and Section 4 of this paper; 4. Patterson et al. 1995; 5. Jones & Watson 1992. 



physics is required. For VW Hyi, WW2002 found that the 
observed rate P = 7.5 x 10~*, which is deduced from an in- 
crease of P from 23 s to 90 s in 24 hours (Warner &: Woudt 
2006; hereafter WW2006), may be explained by the angular 
momentum extracted from the equatorial belt as it feeds into 
the expelled gas. A concomitant precipitate drop in EUVE 
flux during just this phase (Mauche, Mattel & Bateson 2001: 
hereafter MMB; WW2002) is direct evidence for the reduc- 
tion of mass falling onto the primary. 

In Table 2 the large increase of P in the last observed 
DNOs in OY Car is direct evidence for rapid deceleration of 
the belt - but no EUV flux has been observed at this stage. 
Also, a large temporary reduction in EUV flux from SS Cyg 
is seen late in outburst (MMB), and it is during this phase 
of outburst that Jones & Watson (1992) found that a large 
increase in P, rather than a decrease in phase coherence, can 
explain the range of observed DNO phases. The one very late 
phase EUV flux observation in U Gem (MMB; Long et al. 
1996) also shows a precipitate drop, starting midway down 
the optical decline as in VW Hyi, but no DNOs have been 
detected at this time. The increase of P after T = 13 d in 
TY PsA late in the superoutburst of June 1984, although 
not greatly in excess of what is predicted by equation (2), 
is similar to what is seen in OY Car. These results for the 
late declines of VW Hyi, SS Cyg, OY Car, U Gem and TY 
PsA are indications that propellering is a common event in 
late stages of at least some dwarf novae outbursts. 

Although there is not yet a complete theory of accretion 



onto the primary during the late stages of outburst we can 
provide an illustrative example of why the expansion of the 
magnetosphere may overtake the expansion of the corotation 
radius, leading to propellering. Prom the definition of the 

corotation radius we have 

drco/dt = 1.00 X 10*P"^/^PMy^(l)cms"\ (3) 

The radius of the magnetosphere, which we take here to 
be effectively ro, is, for a dipole field, (e.g. Warner 1995) 

ro = 5.15 X 10'VM^Af-'''''(l)M-^/'' cm, (4) 

where is the magnetic moment of the primary (or, in this 
case, the equatorial belt), i.e. n = B{h)B?(l), where B{b) is 
the field in the belt. If we use the theory of field enhancement 
given in equation (14) of WW2002, and adopt appropriate 
values for the parameters, we find 

B3{h) ~ l.ld,BV^{l)Mli\h)p-''''M-^'\l)G, 

which leads to 

ro = 1.71xlO^B^/'^(l)M2Y''(6)P"^/^M-^/''(l)M-^/'^cm, (5) 

which can be differentiated to provide dro/dt. Before pro- 
gressing further, though, we need an estimate of M and in 
particular the way it varies in the late decline stages of out- 
burst. There is as yet no theory for this, but in the case of 
VW Hyi we can use the EUV flux, before propellering starts, 
as a proxy for M. From MMB (see fig. 1 of WW2002) we 
find M oc exp (-t/0.36 d). Adopting Mi(l) = 0.6, Mw ~ 5, 
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M22(6) = 5, Br{l) = 10"* G and P from eq. (2), then 
eq. (5) gives dro/dt ~ 2.2 x lO^cmd"^ and eq. (3) gives 
drco/dt ~ 3 X lO^cmd"^. These rough estimates indicate 
that propellering can indeed be a consequence of the rapid 
drop in near the end of outburst. The propellering will 
continue until the equatorial belt is decelerated to the point 
where equilibrium accretion from the inner edge of the disc 
is re-established. 

As shown by Spruit & Taam (1993; hereafter ST), un- 
less LOs is very large, very little gas is ejected from the sys- 
tem - it is merely moved outward in the disc and will be 
accreted later. In the case of VW Hyi, accretion at a low 
level (Ml 6 ~ 1) is frustrated for about 24 h while it is being 
centrifugally expelled by magnetic coupling to a belt that 
was accumulated during several days at Mie ~ 50; the an- 
gular momentum drain requires only that the rejected gas 
be moved out to ~ 10J?(1). The analysis by ST shows that 
diffusion in the disc brings some of the expelled gas back to 
the inner edge, with the result that the increasing pressure 
there reduces the radius of the inner edge and the condition 
Ws < 1 is eventually regained. In the case of dwarf novae this 
is not a cyclical process (as investigated by ST) because the 
inner disc is rapidly drained of mass after the propellering 
ceases. 



4 COHERENCE OF DNOS 

From the earliest years of their study the DNOs have been 
recognized to have variations in period or phase that occur 
on short times scales, as well as through an outburst. Super- 
imposed on the general spin-up and spin-down arc sudden 
small changes of period. In general the DNOs are more sta- 
ble at maximum of outburst, and deteriorate in stability as 
the outburst progresses - this was first made quantitative by 
Hildobrand, Spillar & Stiening (1981) from optical measure- 
ments and by Cordova et al. (1984) from X-ray observations. 
The disappearance of DNOs (if present at all) from Fourier 
transforms of light curves in the late stages of dwarf nova 
outbursts has been attributed to this efi^ect. 

In the absence of a physical model for DNOs it was 
a natural course to examine their behaviour by statis- 
tical analyses (see Chapter 8 of Warner 1995a for an 
overview) and comparison with particular models. These in- 
cluded noise-excited damped harmonic oscillator (Robinson 
& Nather 1979) and random walk in phase of a sinusoidal 
oscillator (Horne & Gomer 1980; Cordova et al. 1984). These 
were the best that could be achieved for the low amplitude 
oscillations, but examination of the few examples of largo 
amplitude oscillations has repeatedly thrown doubt on the 
value of such an approach. Relevant statements are: (refer- 
ring to phase measurements during an outburst of CN Ori) 
"the remarkable feature of this diagram is the degree of con- 
tinuity in the phase points, indicating that the oscillations 
are coherent throughout the length of the run" (Warner & 
Brickhill 1978); (from analysis of soft X-ray pulsations in 
SS Cyg during outburst) "the random walk in phase model 
does not offer a simplified description of the data" (Jones 
& Watson 1992); "when best defined, the DNOs show the 
presence of an underlying clock with abrupt changes of pe- 
riod, on which the phase noise seems never to cause loss of 
knowledge of phase" (Warner 2004). The set of phase mear 



Table 3. Observing runs on TY PsA 



Run Date (1984) V T (d) Unt (s) Duration (h) 



S3370 25 June 13.2 13 3 3.7 

S3372 26 June 13.9 14 3 2.7 

S3377 27 June 14.9 15 3 0.7 

S3409 22 Sept 13.0 1 4 2.9 

S3412 23 Sept 13.9 2 4 2.0 

S3412 24 Sept 14.9 3 4 2.5 



surements that best displays the basic properties of DNO 
phase variations is that in fig. 10 of Warner & Woudt (2002), 
where it can bo seen that the departures from a smoothly 
varying phase (caused by the general increase of period dur- 
ing decline) arc actually short sections of coherent oscilla- 
tions caused by abrupt but temporary changes of period. We 
shall see below that there is a simple model that explains 
this behaviour. 

To prepare for this we need to describe further phenom- 
ena to those already introduced in Section 2. Additional to 
the classic DNOs there are two kinds of quasi-periodic os- 
cillation (QPO), one kind having periods typically fifteen 
times larger than the DNOs, and another type with peri- 
ods ~1000-3000 s (see table 1 of Warner 2004). On rare 
occasions double DNOs are observed, with their difference 
frequency equal to the QPO frequency. The LIMA model 
explains classic DNOs as reprocessing of the beam of high 
energy radiation (emerging from the accretion zone) from 
the accretion disc, and the additional DNO as reprocess- 
ing from a prograde traveling wave in the inner disc which 
is itself responsible (via obscuration and/or reprocessing of 
radiation from the central source) for the QPOs seen in the 
light curve. We therefore refer to the direct DNO as having a 
sidereal period, and that reprocessed from the traveling wave 
as having a synodic period (WW2006). On occasion VW 
Hyi shows alternation between the two types of oscillation 
(fig. 12 of WW2006), but has other variations (frequency 
doubling and tripling) that complicate the situation, so we 
here investigate a simpler system. 

4.1 Reanalysis of the light curves of TY PsA 

The dwarf nova TY PsA was observed during a superout- 
burst in June 1984 and showed DNOs of sufficient amplitude 
for them to be seen directly in the light curve. These gave 
the finest opportunity to that time for following the DNO 
phase variations on individual nights, and lent support to 
the notion of decreasing coherence through outburst decline 
(Warner, O'Donoghue & Wargau 1989; hereafter WOW). 
Here, in the light of later experience, we reinvestigate some 
of the DNOs present in that superoutburst and also the 
DNOs present in a normal outburst three months later. The 
relevant runs are listed in Table 3. More complete details 
and the light curves can be found in WOW. 

It will be seen that the brightness of TY PsA fell at 
almost exactly the same rate at the ends of the two out- 
bursts; this is a well known feature of dwarf novae - normal 
and sup erout bursts have nearly identical decline light curves 
(Warner 1995a: and see WW2002 for details of VW Hyi light 
curves). 
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Figure 1. The phase diagram of run S3370 of TY PsA. Numbers 
along the top of the diagram are the DNO periods (in seconds) 
measured for the sections of the light curve indicated by dotted 
vertical lines. The time ranges of these sections were defined so 
that times when O — C is systematically increasing are isolated. 
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Figure 2. The O — C diagram of a ~ Ih section of run S3372. 
O — C can be seen to alternate between increasing and decreasing. 



coherence and hence invisibihty in the FT. However, using 
an j4/c1? search, with a DNO period initially extrapolated 
from the P trend on previous nights, has revealed that there 
is one section of the light curve, lasting for almost exactly an 
hour, that does have a detectable signal. This is shown as an 
A/$ diagram in Fig. 2, with reference period 38.0 s, where 
it is seen to be another example of alternation between two 
distinct periods, in this case ~ 36.96 s and ~ 39.93 s, which 
have a beat period ~ 500 s, similar to the previous night, 
but now the switching time scale has doubled to ~ 1000 s. 



We have analysed the light curves in the standard 
fashion using Fourier transforms (FTs), and also ampli- 
tude/phase (j4/$) diagrams, where we fit sinusoids, with 
periods derived from the FTs, by least squares to short sec- 
tions of the light curve to determine the amplitudes and 
phases of the DNOs as functions of time. 

4.1.1 Run S3370 

In WOW the A/^ diagram was (correctly) described as hav- 
ing coherent oscillations of only short duration, responsi- 
ble for the noticeable broadening of the peak in the FT, 

compared with DNOs of longer coherence in three previous 
nights. As plotted in WOW the A/^ diagram appears quite 
chaotic, but we noticed that with suitable adjustments of 
sections of the diagram, permitted by the 2ti phase ambi- 
guity of any point, a more systematic picture appears. We 
show this in Fig. 1, where $ is measured relative to a fixed 
period of 30.17 s, and where we have partitioned the dia- 
gram into regions of clearly increasing or decreasing $ and 
have indicated the DNO periods P derived from FTs of each 
of those short sections. It is clear that P alternates between 
two values - approximately 30 s and 31 s, and it does so 
on a timescale ~ 0.011 d, which is close to the beat period 
950 s between the alternating values. An FT of Fig. 1 
shows a strong peak at ^ 500 s, with a sub-harmonic at 
^ 1000 s. Much of the $ variation in Fig. 1 is therefore de- 
scribable by a switching between a sidereal period ~ 30 s 
and a synodic period ~ 31 s, and back again, at the beat 
period between the two. The similarity to the behaviour of 
VW Hyi, mentioned above, is evident. 

4.1.2 Run S3372 

In WOW it was stated that Run S3372 does not contain 
any DNO signal, which, if it were the case that DNOs were 
really present, was attributed to a probable complete loss of 



4.1.3 Run S3377 

WOW did not detect any DNOs in this short run, and it is 
true that an FT of the whole run docs not show any convinc- 
ing signal, but more detailed investigation shows an initial 
short-lived presence of quite high amplitude but uncertain 
period, and a persistent DNO at 47.0 s of low amplitude in 
the second half of the run. At this stage TY PsA was still 
about 1 mag above its quiescent magnitude. 



4.1.4 Run S3409 

In WOW it was noted that this light curve possesses a DNO 
of large amplitude with a period of 26.5 s. Its properties are 
similar to those seen for run S3365 in fig. 13 of WOW and 
are characteristic of the early outburst decline, rather than 
late decline, and are not of particular interest here. 



4.1.5 Run S3412 

WOW noticed only a QPO with period ~ 245 s in this run. 

However, we find that there is also a DNO present that 
changes its period so rapidly that its peak in the FT of the 
raw light curve is of low amplitude and greatly broadened, 
as seen in the FT displayed in the left hand panel of Fig. 3. 
The DNO is clear in the A/^ diagram (right hand panels of 
Fig. 3), where the reference period is 38.59 s, but we have 
not found it possible to find an unambiguous interpretation 
of the variations. What is evident, though, are at least four 
large peaks of amplitude, spaced roughly 500 s apart and 
plotted as larger symbols in Fig. 3, and concomitant rapid 
switches between different periods of roughly 38.6 s and a 
period probably about 2 s longer. We have the impression 
that while the ~ 500 s modulation (at twice the QPO period) 
is present the behaviour is similar to what is seen in Fig. 2. 
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Figure 3. The FT (left hand panel) and O — C and amplitude curves (right hand panels) of run S3412 of TY PsA. O — C is calculated 
relative to a fixed period of 38.59 s. The DNO shows up in the FT as excess power between roughly 0.023 and 0.030 Hz. Large symbols 
are used to emphasise four pulses in the amplitude curve. 



4.1.6 Run 33417 

Run 3417, unremarked by WOW, has a QPO with period 
~ 241 s, which is similar to that seen on the previous night. 
We do not detect any DNOs in this run. 

4.1.1 Discussion 

We have found that DNOs are present later in the outburst 
light curve of TY PsA than had previously been realized. 
There is abundant evidence of the kind of switching between 
sidereal and synodic DNOs only previously seen in VW Hyi 
(WW2006). On the outburst dechne a QPO periodicity at 
~ 250 s and its multiples is present. From a different per- 
spective, we can say that a period ~ 1000 s and its first and 
second harmonics are favoured, which again resembles the 
behaviour of VW Hyi (WW2006). 

In the hght curves of TY PsA analysed here we have not 
found any evidence of IpDNOs, having typically four times 
the periods of DNOs (Warner, Woudt & Pretorius (2003; 
hereafter WWP). However, they are occasionally present in 
this star, as shown in fig. 3 of Warner (2004). 



5 A POSSIBLE CONNECTION BETWEEN 
THE MULTIPLE DNOS IN WZ SGE AND 
LONGER PERIOD DWARF NOVA 
OSCILLATIONS 

We aim in this Section to ofTer a model for the site of the 
IpDNOs alternative to that put forward by WWP. Our ap- 
proach is via a hitherto unexplained periodic modulation in 
the dwarf nova WZ Sge. 

5.1 Observational overview 

WZ Sge is unusual in possessing a set of DNOs even during 
quiescence, which show complicated behaviour during the 
very rare outbursts. The key phenomena are: 



Quiescence 

(i) Two optical modulations are observed, at 27.868 and 
28.952 s, discovered by Robinson, Nather & Patterson 
(1978). In general either one or the other modulation is 
present, but occasionally both are seen simultaneously. 

(ii) No harmonics of the 27.87-s modulation are observed, 
but on rare occasions a first harmonic of the 28.95-s modu- 
lation is seen (Provencal & Nather 1997). 

(iii) The 27.87-s signal has substantial phase noise but 
is a good long-term clock (Patterson 1980; Patterson et al. 
1998). The long-term stability of the 28.95-s signal has not 
been established. 

(iv) Hard X-ray, soft X-ray and UV modulations have 
been observed at the 27.87-s period (Patterson et al. 1998; 
Skidmore et al. 1999); however, the 28.95-s periodicity is not 
associated with any such high energy radiation. 

(v) Sidebands are observed to the 27.87-s modulation, 
which have been used in support of a rotating magnetic 
model for the white dwarf in WZ Sge (Warner, Tout & 
Livio 1996; Lasota, Kuulkers & Charles 1999; Patterson et 
al. 1998; WW2002). 

(vi) Patterson (1980) suggested that the 28.95-s signal 
could be generated from the 27.87-s signal if the latter is 
reprocessed from a site moving revolving progradely in orbit 
at the beat period between the two, i.e. with a period of 
742 s. W2002 found just such a period in light curvei{3 of 
WZ Sge and associated it with the same type of traveling 
wave that generates QPOs in the LIMA model. 

Outburst 

(i) DNOs were observed during the 2001 superoutburst 
of WZ Sge, with a dominant period of ~ 15 s and short- 
lived modulations at ~ 6.5 s thirty days after maximum 
(Knigge et al. 2002), and ~ 18 s fifty days after supermax- 
imum (Welsh et al. 2003). The 27.87-s periodicity had not 

1 The statement by Welsh et al. (2003) that the 742-s periodicity 
is very rarely seen is not correct: the dips that it generates in the 
light curves of WZ Sge are very rarely absent. 
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reappeared two years later (Mukai & Patterson 2004) . It was 
undetectable for nearly 18 years after the 1978 outburst. 

(ii) In contrast to the latter, the 28.95-s modulation was 
detected by these observers with an apparently unchanged 
period at various times during the outburst. 

These behaviours show certain similarities to other 
dwarf novae during outburst - namely, the increase in P 
on the dechne of outburst and the frequency doubling ob- 
served at high M in SS Cyg (Mauche 1998), which may 
be caused by changing geometry (WW2002) or by the in- 
creasing importance of higher magnetic multipoles close to 
the primary. Furthermore, the independence of the 28.95-s 
period on luminosity of the system resembles that observed 
for the IpDNOs, which have been interpreted as arising from 
accretion directly onto the primary, i.e., associated with the 
spin of the primary, not of its equatorial belt (WWP). 

There are also interesting comparisons to be made with 
GK Per, an IP nova remnant showing dwarf nova outbursts. 
It has a largo amplitude 351. 3-s modulation in the X-ray re- 
gion during outburst (Watson, King & Osborne 1985; Hel- 
lier, Harmer & Beardmore 2004) which persists during qui- 
escence (Norton, Watson & King 1988). The 351-s period 
is detected in optical light curves in quiescence (Patterson 
1991) but not convincingly during outburst (cf. Nogami et 
al. 2002). Optical quasi-periodic oscillations are seen in qui- 
escence and outburst with a mean period near 380 s (Pat- 
terson 1981, 1991; Mazeh et al. 1985; Nogami et al. 2002). 
A longer period QPO is seen during outburst near 5000 s, 
which is the beat period between the 351-s and 380-s mod- 
ulations, in X-rays (Watson et al. 1985; Hellier et al. 2004) 
and optical emission line profiles (Morales-Rueda, Still & 
Roche 1999) but not during quiescence. The 5000-s X-ray 
QPOs have been explained by Hellier et al. (2004) in terms 
of the LIMA model, with a traveling wave at the inner edge 
of the disc. 

During outbursts of GK Per there are no observed 
DNOs (which would be expected to progress to lower peri- 
ods than 351 s if they existed). Complete absence of DNOs 
is a property of IPs in general and is interpreted as implying 
that the magnetic field on the primary is strong enough to 
prevent slippage of the equatorial belt (WW2002). 

Finally, we note that WZ Sge has probable low am- 
plitude 27.87-s modulation in soft X-rays and even weaker 
in haxd X-rays (Patterson et al. 1998), the DNOs in dwarf 
novae and nova-like systems have large amplitude soft X- 
ray modulation, but at most weak signals in hard X-rays, 
and GK Per, as with all IPs, has predominantly hard X-ray 
modulation. 



5.2 Generation of the 28.95-s modulation in WZ 
Sge 

The observational properties listed in the previous section 
show that WZ Sge appears to be a bridge between the stan- 
dard IP and the LIMA models. In quiescence the primary's 
magnetic field is barely strong enough to couple the accret- 
ing surface regions to the interior, producing considerable 
phase noise in the 27.87-s spin modulation, but during out- 
burst a slipping accretion zone is possible. Two possible sites 
suggest themselves for the generation of the 28.95-s peri- 
odicity: (a) the main body of the primary, rotating with 



that period, in which case 27.87 s must be the period of a 
permanently slipping equatorial belt, or (b) the region just 
outside the corotation radius in the disc. Both sites retain 
their rotation periods independent of the very large mass 
transfer variations that occur during outburst. The former 
is unlikely because we would expect X-ray production from 
accretion onto the primary as well as its equatorial belt, but 
no 28.95-s X-ray modulation is observed. The latter holds 
more promise. 

For WZ Sge, 28.95 s/27.87 s = 1.039. At this point we 
mention another rapid rotator, AE Aqr, which has P(l) = 
33.078 s (Patterson 1979) and such a low M that almost all 
accretion is propollorod out of the system (e.g. Eracleous 
& Horne 1996). In his discovery paper Patterson (1979) 
pointed out that, as well as the prominent 33.08-s modu- 
lation, there were frequent appearances of a ~ 34.6-s modu- 
lation and its first harmonic. The more comprehensive later 
study (Patterson 1994) shows that the harmonic of this fea- 
ture, with P ~ 17.4 s, is present strongly during 'flaring' 
episodes of AE Aqr, and is only very weakly present in qui- 
escence. Skidmore et al. (2003) also noticed a longer period 
(34.23 s) during a flare. The ratio 34.8s/33.08s = 1.052 is 
comparable to what is observed in WZ Sge and again may 
imply that activity in the vicinity of the corotation radius 
is generating a modulation. In the case of AE Aqr it seems 
certain that propellering, which begins to operate at and 
beyond fco, is a dominating mechanism, so we suspect that 
milder propellering is in action in WZ Sge, contrifuging gas 
out on held lines that connect to the disc just beyond rco, 
during both quiescence and outburst. 

Propellering has already been invoked in WZ Sge, to 
explain the very long intervals between superoutbursts a 
result of the inner parts of the disc being cleared of gas by the 
primarys magnetosphere, preventing outbursts from starting 
near the inner edge of the disc, and allowing large amounts 
of mass to accumulate in the outer regions (Matthews et al. 
2007). 

We are now ready to propose a specific model for the 
generation of modulations from a region near the corotation 
radius. Magnetic field lines connecting from the primary to 
gas revolving at Koplorian velocities in an annulus just be- 
yond Tco are twisted back and produce a retarding torque 
on the primary. In the annulus a blob of threaded gas is con- 
comitantly accelerated azimuthally and radially outwards by 
the field lines, the maximum effect occurring where the field 
(for a tilted rotating dipole) is largest, producing a radial 
velocity Vrei ~ 47/3cs, where 7 is the ratio of the toroidal 
to the vertical field strengths (close to zero at rco), Cs is 
the thermal sound velocity and f) is of the order of unity 
(Aly & Kuijpers 1990, hereafter AK). Immediately outside 
rco the radial velocity is subsonic but twisting of the field is 
already underway. Eventually the threaded gas reaches su- 
personic velocity and shock waves begin to form. AK show 
that the short-lived shocks spiral outwards and remove ex- 
cess angular momentum. For a tilted dipole field a group 
of such shocked areas will develop each side of the disc, be- 
ing carried around at a period close to the Keplerian period 
of the annulus, not far outside of r^o, in which the shocks 
occur. Heating by the shocks will thicken the disc locally; 
it is the cessation of this heating when the field lines snap 
and reconnect that was hypothesized by WW2002 to be the 
excitation mechanism for generating the slow prograde trav- 
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eling wave at the inner edge of the disc, which produces the 
QPOs. 

We suggest that the two shock regions, revolving as 
bright regions in the disc (which individually may be viewed 
differently, depending on orbital inclination, optical thick- 
ness of the disc, etc.) generate the fundamental and/or har- 
monic of the 28.95-s and 30.5-s modulations seen respec- 
tively in WZ Sge and AE Aqr. 

An estimate of the total luminosity Lgh expected from 
the shock regions is provided by AK. Their equation (25) 
can be written 

Lsh = 6.3 X 10"^^ fiiiT-g'' HgCfj crgs^^ (6) 

where rg — rco/10"cm, Hs is the scale height of the disc 
in units of 10* cm, and = Cs/(10^ cms^^). From mod- 
eUing the outbursts of WZ Sge, Matthews ct al. (2007) find 
they require 0.5 < /X31 ^ 3. For Afl(l) = 1.2 (Spruit 
& Ruttcn 1998) wc have rg — 1.50, and therefore Lsh ^ 
1 X lO^^ergs"^. When compared with the accretion lumi- 
nosity of WZ Sge, L ~ 3 X 10*° erg s"^ (Smak 1993) it is 
evident that a detectable weak modulated signal could well 
be generated from the shock regions. The period Psh of this 
signal will not be completely stable - variations of condi- 
tions in the threading annulus will result in small changes 
in Vrei and hence also in Psh- 

We note that according to equation (6) Lsh may ex- 
pected to be observable only in systems with low accretion 
luminosity, as is the case for WZ Sge (M ~ 2 x 10"'''' gs""'; 
Smak 1993) and AE Aqr (M - 1 x 10"gs"^ onto the pri- 
mary; Wynn, King & Horne 1997). At higher M, H only 
increases oc M^/^° and so Lsh /Lace could be expected to 
fall rapidly during the outburst of a dwarf nova (or in the 
high state of a nova- like variable). But the possibility of an 
enhanced field in the equatorial belt, leading to an increase 
of by a factor of a few (WW2002) might, from equation 
(6), still give rise to a detectable Psh signal. In that case, it 
is possible that all the IpDNO signals that arc observed in 
CVs arise not from accretion zones on the body of the pri- 
mary, as suggested in WWP, but from the shock regions just 
outside the corotation radius. These two sites may be obser- 
vationally distinguishable by photometry of IpDNOs during 
eclipse. The nova-like EC 21178-5417 could provide this op- 
portunity - it has IpDNOs with periods ~ 95 s and deep 
eclipses (WWP). For Mi(l) = 0.8 the source of the modular 
tion is either at -R(l) 7 x 10^ cm or near rco ~ 3 x 10® cm, 
a difference that might be disccrnable during eclipse ingress 
or egress, though it must be kept in mind that the whole disc 
could be involved in the optical signals through reprocessing 
of radiation from the sources. 



6 CONCLUSIONS 

A principal result of this study is further evidence that the 
variations in period or phase of DNOs are best described by 
switching between two periods - the sidereal period, deter- 
mined by the rotation period of an accreted equatorial bolt, 
and a synodic period, probably generated by reprocessing of 
the radiation emitted from accretion zones on the equato- 
rial belt. This effect is only seen where the DNO amplitudes 
are large enough to enable detailed analysis, but may be 
suspected to hold in general. 



We point out that propelloring of accreting gas may 
commonly occur towards the end of dwarf nova outbursts, 
depending on the existence of a sufficiently strong magnetic 
field on the primary. We suggest that the longer period 
DNOs might be located at regions just beyond the corota- 
tion radius, rather than on the surface of the primary, and 
that previously observed but not understood periodicities in 
WZ Sge and AE Aqr arise from this source. 
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